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O ■ Abstract 

By considering the potential parameter F as a function of another potential parameter 
A [47], We successfully extend the analysis of two-dimensional autonomous dynamical 
^1^' system of quintessence scalar field model to the analysis of three-dimension, which makes 

^ ' us be able to research the critical points of a large number of potentials beyond the 

exponential potential exactly. We find that there are ten critical points in all, three points 
' P'i,b,& are general points which are possessed by all quintessence models regardless of the 

^r~i ' form of potentials and the rest points are closely connected to the concrete potentials. It 

' is quite surprising that, apart from the exponential potential, there are a large number of 

\ potentials which can give the scaling solution when the function /(A)(= r(A) — 1) equals 

' zero for one or some values of A* and if the parameter A* also satisfies the condition 

00 ! Eq.(16) or Eq.(17) at the same time. We give the differential equations to derive these 

potentials V{(j)) from /(A). We also find that, if some conditions are satisfied, the de- 
Sitter-like dominant point P4 and the scaling solution point or Pio) can be stable 
^ \ simultaneously but Pg and Pio can not be stable simultaneity. Although we survey scaling 

5^ \ solutions beyond the exponential potential for ordinary quintessence models in standard 

general relativity, this method can be applied to other extensively scaling solution models 
studied in literature [46] including coupled quintessence, (coupled-)phantom scalar field, 
k-essence and even beyond the general relativity case oc prp. we also discuss the 
disadvantage of our approach. 
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1 Introduction 



Scalar fields had played an essential role in modern cosmology in the past semi-century. This 
assumed scalar field had been used for various purposes in different cosmological research 
aspects[l], such as to drive inflation, to explain a time variable cosmological "constant" and 
so on. Especially, after the discovery of the accelerating expansion of universe, it has gained 
another hotly discussion as the candidate for dark energy. There are so many scalar fleld dark 
energy models, such as quintessence model[2-13], non-canonical scalar fleld model ( including 
K-essence[14-17], phantom [18-22], B-I scalar fleld[23-29] and so on) and coupled scalar fleld 
model[30-31]. There are also detailed studies on the multi-scalar fleld models which give 
an effective state equation Weff passing through the phantom divide line {w = —1) [32-34]. 
Some of these multi-scalar fleld models [35-36] can always evolve onto the regime of scalar fleld 
dominance Ag//^ > 37 even if each fleld has too steep a potential to drive the accelerating 
expansion. For all of these scalar fleld models we mention above, the important thing is to 
choose different form of kinetic terms and different potentials from a fundamental physical 
motivation or directly from the observation. As are expected, these different scalar fleld 
models will give different cosmological evolutions, different evolutions of state equation w, 
different values of sound speed and different cosmological perturbation. So they can in 
principle be distinguished or excluded by the increasing observation data. 

The phase-plane analysis of the cosmological autonomous system is an effective method 
to flnd the cosmological scaling and dominant attractor solutions. A phase-plane analysis 
of cosmologies containing a barotropic fluid and a scalar fleld with an exponential potential 
was presented [3 7]. Hao and Li studied the attractor solution of phantom scalar fleld with 
the exponential potential [38-39]. On the other hand, L.Amendola considered the case of 
coupled quintessence[31]. The case of phantom scalar fleld interacting with dark matter 
was also investigated [40-41]. Guo also investigated the properties of the critical points of 
multi-field model with an exponential potential [42-43] One may realize that the potentials 
investigated in all these papers are the exponential form. Disregarding the important roles 
of the exponential potential in higher-order or higher- dimensional gravity theories and string 
or kaluza-klein type models, the reason that why they are choosing the exponential potential 
may be that, only the exponential potential can give a two-dimension autonomous system. 
Since in this case the value of the parameter F equals 1 and then another parameter A equals a 
constant (see Eq.(4) for the definition of parameters F, A ), so the system(see Eqs.(5-7) below) 
will reduce to the two-dimension autonomous system. However, authors also considered the 
more complicated case when A is a dynamically changing quantity[44-46]. They applied 
the discussion of constant A to this case and obtained the so-called "instantaneous" critical 
points. For example, if F is a constant (but does not equal one), saying F = (n + l)/n, 
the corresponding potential is the inverse power-law potential V{(f)) = Vo(j)~"' with n > 0. 
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One of the critical points {xc,yc) = {^/VQ, [1 — will become the "instantaneous" 

critical point {x{N) = \{N)/V6, y{N) = [1 - A(iV)V6]i/2)_ ^^en T > 0, A(A^) will decrease 
toward zero, then the "instantaneous" critical points will eventually approach x{N) and 
y{N) — s> 1. This method is not exact here and obviously the critical point is not a true critical 
point. Recently, a solution of multiple- attractor in three- dimension autonomous system of 
the quintessential models was studied in literature [47]. After writing the parameter F as a 
function of A, the author obtained a tracker solution which is different from those discovered 
before and found a solution of multiple-attractor. Here we will extend the idea to an arbitrary 
function r(A). We will find out all the critical points of the dynamical autonomous system, 
and then investigate the properties of the critical points and their cosmological implications 
in general. Regarding parameter F as a function of A is a quiet efficient approach since 
we can investigate many quintessence models with different potentials. Giving a concrete 
form of function r(A) is equivalent to give a concrete form of potential V{(j)) since we can in 
principle figure out the potential via the relation between parameter F and A. What are the 
general properties of the critical points when we consider the three-dimension autonomous 
system? Does there also exists scaling solution when we consider any function of F(A)? 
Among all the critical points which critical points are the critical points for all quintessence 
and which are only relative to the concrete potentials? In our paper, we will try to shed light 
on these issues. The paper is organized as follows: in Section 2 we present the theoretical 
framework and give the differential relation between the function F(A) and potential V{(j)). 
We find out all the critical points and investigate their properties in Section 3. We try to 
give the cosmological implications of these critical points in section 4. We briefiy display our 
conclusions in section 5. 



2 Basic theoretical frame 

We start with a spatially fiat Friedman-Robertson- Walker universe containing a scalar field 
(j) and a barotropic fiuid (with state equation pf, = WbPb)- To simply, we give the Einstein 
equations directly: 

H' = ^[\<P' + V{<P)+Pb] (1) 

^ = -Y[0' + (1+w;6)P6] (2) 
The motion equation of the scalar field (p is: 

dV 

4> + 3H4> + — = (3) 
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Following [48], we define the following dimensionless variables: 

Where V = dV{(f)) /d(f),V" = £V {(p) / dcp'^ . Using Eq.(4), Eqs.(l-3) can be rewritten in the 
following dynamical form[37, 46, 48]: 

1^ = -3x + ^Xy' + - w,)x' + (1 + w,){l - y')] (5) 

^ = -^Axy + ^y[{l - w,)x' + (1 + ^,)(1 - y')] (6) 

^ = -V6X\r - l)x (7) 

where N = ln{a). Here we should emphasize that Eqs.(5-7) is not a dynamical autonomous 
system since the parameter F is unknown. However, if we consider F as a function of A, 
namely 

F(A) = /(A) + 1 (8) 

then Eq.(7) becomes: 

^ = -V6AV(A)x (9) 

Hereafter, Eqs(5-6) and Eq.(9) are definitely a dynamical autonomous system. We will see 
that F as a function of A can cover many quintessential potentials. The three-dimension 
autonomous system reduces to two-dimension autonomous systems when /(A) = 0. In 
this case, the potential is the exponential form which has been completely studied in many 
literatures. When /(A) equals a nonvanishing constant fx, then the potential is proportional 
to (ci0 -|- C2)~^^-^^, which is just the potential which has been considered as "instantaneous" 
critical points[48]. Generally speaking, we can analyze any explicit function. For some 
more comphcated form, F(A) = 1 + ^ - ^ corresponds to \/(</)) = i^^^^^^, F(A) = 1 ± ^ 
corresponds to V{(f)) = Voe^''^^^+f^'>/^ , F(A) = 1 + ;^ corresponds to V{(f)) = V^e^l'^. The form 
of F(A) = 1 -|- -|- y, which corresponds to = (^pp^s^p-, was considered as an interesting 
cosmological model where the universe can evolve from a scaling attractor to a de-Sitter-like 
attractor by introducing a possible mechanism of spontaneous symmetry breaking[47]. 

In the paper [47], the author gave an approach to obtain the potential ^(0) as follows: 
Since the potential V^(0) is only a function of the field 0, then the parameters A and F can 
be written as a function of field: A = -P(0), F = Q{4'). If the inverse function of P{4>) exists, 
then we have: 

T = Q{P~\X)) ^ nX) (10) 
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Using the definition of A and T, V" can be written as V" = = F{V,V'). Let 

h = V', then 

Now Eq.(ll) is a one-order differential equation of h and V. Figuring out h{V), the potential 
can be solved from equation V'{(f)) = h(y{(j))). 

Here we introduce another easier approach to get the potential V{(j)). We start with 
§ = ^W = -^^^177 = -yr^^^^- Using the definition of A and T, and the Eq.(8), 
we get a one-order differential equation of A and V: 

Integrating out A = \{V), using the definition of A, then we have following differential 
equation of potential: 

-d<j) (13) 



VX{V) 

So Eq.(12) and Eq.(13) give the route to obtain the potential V = V{(j)). As far as we 
know, there are too many investigations to the two-dimension autonomous system (where 
the potential is exponential) but have not general investigations to the dynamical properties 
of three dimensional autonomous system. It is maybe very interesting to consider this issue. 
We know that previous exactly analysis to the critical points of the quintessence model are 
based on a concrete form of potential (i.e., the exponential form). In this case it is not easy 
to distinguish which critical points are common to all the quintessence models and which 
are only related to the special potentials. In view of what we mention above, we will take a 
new route in next section to investigate the critical points of the autonomous system with 
an arbitrary function of /(A). Furthermore, the results can be easily applied to any other 
concrete potentials as long as they can be solved from Eqs. (12-13). 



3 Critical Points and their Properties 

It is easily seen from Eq.(9) that A = 0,x = or /(A) = can make dX/dN = respec- 
tively. The critical points listed in TABLE 1 can be found from the Eqs. (5,6,9) after setting 
dx/dN = dy/dN = dX/dN = 0. The properties of each critical point are determined by the 
eigenvalues of the Jacobi matrix of the three-dimension autonomous system. For a general 
three-dimension autonomous system: 

^ = fi{x,y,X) 

^ = f2{x,y,X) (14) 
j^ = f3{x,y,X) 
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The function /i, /2 and /a are only the function of x, y, A, no variable N and other 
variables, we call this dynamical system as autonomous system. If /i, /2 and /s are only 
a linear combination of x,y,X, Eq.(14) is linear autonomous system. Its critical points 
(xc, He, Ac) can be found from the set of functions /i = /2 = /s = 0. Obviously, Eqs. (5,6,9) is 
not a linear autonomous system. However, the local behavior of the nonlinear autonomous 
system near a critical point can be deduced by linearizing the nonlinear system about this 
point and be studied using the linear autonomous system analysis method. The properties 
of each critical point are determined by the eigenvalues of the Jacobi matrix A, where 



For a hyperbolic critical point^, if all the eigenvalues of A or the real part of these eigenval- 
ues are negative, the critical point is stable. This is to say, as long as one of the eigenvalues 
or the real part of these eigenvalues is positive, the critical point must be unstable. However, 
if the critical point of nonlinear autonomous system is a nonhyperbolic point^ and the rest of 
its eigenvalues having negative real part, the properties of this point can not be simply de- 
termined by linearization method and need to resort to other more complicated methods[50]. 
From TABLE 1, we can see that point P4 is just this kind of point. In previous literatures [30, 
31, 37, 51], the authors generally neglected this nonhyperbolic point when they met it. In 
fact this point also has the important cosmological implication as other critical points and 
should not be ignored. We will explore the properties of this nonhyperbolic point P4 in 
our paper using the center manifold theorem [50] (The full analysis process is given in the 
Appendix). We list all the points and their properties in the following TABLE 1. Note that 
we have neglected the cases with y < since the system is symmetric under the reflection 
(A, X, y) — > (A, X, —y) and time reversal t —t. 

^ Actually "critical point" in this paper is also called the "equilibrium point" in mathematics or "fixed 
point" in some physical literatures, a hyperbolic critical(equilibrium) point is the critical(equilibrium) point 
which has no eigenvalues with zero real part. 

^ i.e., its eigenvalues exist zero value or have zero real parts. 



A 



dfi{x,y,X)/dx dfi{x,xy,X)/dy dfi{x,y, X)/dX 
df2{x,y,X)/dx df2{x,y,X)/dy df2{x,y, X)/dX 
df3{x,y,X)/dx df3{x,y,X)/dy df3{x,y, X)/dX 



(15) 
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f Ar-. Xr-. Vr^) 


eigenvalues 


Stability 


Pi 


(0,1,0) 


3(1 - Wh),3,0 


unstable node 


P2 


(0, -1,0) 


3(1 _ m),3,0 


unstable node 




(0,0,0) 


-3(1 - m) 72,37/2,0 


saddle point 




(0,0, 1) 


-3, -37,0 


stable node for 
/(O) > 




(Aa,0,0) 


-3(1 -if;fe)/2, 37/2,0 


saddle point 




(A., 0,0) 


-3(1 -tf;fe)/2, 37/2,0 


saddle point 


P7 


(A*, 1,0) 


-V6A2d/„3(l - ^/;fe), i(6 - VeA,) 


saddle point 


Ps 


(A.,-1,0) 


v^A2rf/„3(l-w;6),i(v^A, + 6) 


saddle point 


P9 


(A*, -^A*, -y/l - |A2) 


|(A^-6),A2-37,-A2A.C 


Eq.(16) 


Pio 


/■V v^7 V67(l-«)6)x 


-SX^jdf^, l{wb - 1) 
^ VOp) ^24^ - (97 - 2)Xl 


Eq.(17) 



Where /(O) is the value of function /(A) at A = 0, c//* = ^^^|a*- We limit the range 
of Wfe(= 7 — 1) as < < 1, = for matter and 1/3 for radiation. Aa means an 
arbitrary value and A,,, is the value which makes /(A,,,) = 0. So points Pj^io appear only 
if the function /(A) can be zero for one or more values of A^,. Here we simply consider that 
only one value A* makes the function /(A) zero. 

However, readers should keep in mind that, to make dX/dN = in Eq.(9), we let A = 
0, X = and /(A) = separately, and then find out all the points listed in TABLE 1. But we 
do not consider one special case that A^/(A) 7^ and then dX/dN ^ when A = 0. In this 



case. Pi and P2 are no more critical points. For example, the product A^/(A) 



7^0 



even if A = for the potential V{(j)) = Vo[cosh(a0) — 1] + A. So the necessary condition for 
the existence of equilibrium points with x 7^ is A^/(A) = 0. 

A^ < 6 is the condition for critical point Pg to exist and Eq.(16) is the condition for Pg 
to be a stable node. 

A^ < 37 and X^df^ > (16) 
A^ > 37 is the condition for critical point Pig to exist and Eq.(17) is its stable condition. 



247^/(97 — 2) > A^ > 37 and X^df\ > for Pig being a stable node 
A^ > 247^/(97 — 2) and X^df\ > for Pig being a stable spiral 

the density parameter of (f) field and its equation of state are, respectively: 



(17) 



W0 



(18) 
(19) 
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In order to investigate the expansive behavior of scale factor a, we also represent the decel- 
erating factor: 

(20) 

= ~ Wb)x^ - {1 + Wb)y^ + {wb + I)] 
We list the other properties of these critical points in TABLE 2. 











decelerating factor(q) 


Pi 


(0,1,0) 


1 


1 


2 


P2 


(0,-1,0) 


1 


1 


2 


P3 


(0,0,0) 


Undefined 





(3^6 + l)/2 


P4 


(0,0,1) 


-1 


1 


-1 


n 


(Aa,0,0) 


Undefined 





(Swb + l)/2 


Pe 


(A., 0,0) 


Undefined 





{3wb + l)/2 


P7 


(A*, 1,0) 


1 


1 


2 


Ps 


(A*, -1,0) 


1 


1 


2 


P9 


(A*, ^A=^, -y/l — |A2) 


A^/3 - 1 


1 


A^/2 - 1 


Pio 


(\ Ve-y V67(l-«i6)x 


Wb 


37/A^ 


{3wb + l)/2 



4 Cosmological Implications 

After giving all the critical points and their properties of the three-dimension autonomous 
system, we will investigate their cosmological implications. We will show some interesting 
results which have not been found previously in other literatures. Moreover, we will also 
response to the questions we have proposed in Section 1. Investigating three-dimension au- 
tonomous system instead of the two-dimension autonomous system can help us consider more 
potentials which can not be investigated via two-dimension autonomous system. Moreover, 
from the view of three- dimension system, we can gain a more deeply understanding than 
from the two-dimension system. For example, we will point out which critical points are the 
critical points for all quintessence and which are only relative to the concrete potentials. We 
can find from TABLE 1 and TABLE 2 that: Though the stability of Points Pi^2 does not 
depend on the form of concrete potentials. Points Pi^2 only exist when A^/(A) = at A = 
. Points ^3,5,6 always exist for all quintessence models and their stability are regardless of 
the form of concrete potentials. Point P4 is also the critical point for all quintessence, but 
its stability depends on the form of concrete potentials. Points Pj^io and their properties 
are closely connected to the concrete potentials since the value of A* is determined by the 
form of /(A), points P7-10 are even inexistence if /(A) 7^ for any A. 
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Of all the points, only Points ^3,5,6 are independent of the function /(A). In fact, they 
have the same properties and can be considered as one point. They are saddle points which 
tell us that the barotropic fluid dominated solution (Ac = 0, Xc = 0, yc = 0) where = is 
unstable. However, even though they are unstable, the phase space trajectories may evolve 
in the vicinity of the barotropic fluid dominated solution for a quite long time and then 
leaves this state to approach to the possible future attractor. However, if 7 = 0, these points 
are found to be a stable attractor and can be used to alleviate the relic density problem in 
inflation model[37]. 

Four of the critical points (Pi,2(Ac = 0,Xc = ±1, Uc = 0) and Pj^siK = K, Xc = ±l,yc = 
0)) are all unstable nodes, which correspond to the solutions where the universe is dominated 
by the kinetic energy of the scalar field {Q^ = 1) with a stiff equation of state {w^ = 1). 

In fact, we can conclude above results with one brief sentence (see TABLE 1): all the 
critical points with yc being zero are not stable points. It tells us that, under the potential 
we considered here, the cosmological solution with the potential energy eventually evolving 
to zero will never be the final state of our universe. This is a quite interesting result since 
we know that the universe will never undergo a regime of accelerating expansion if there is 
no potential energy in quintessence models. 

Therefore, there are only three critical points ^4,9, 10 which correspond to possible late- 
time attractor solutions. We will study their properties and cosmological implications in 
more detail. 

Points P^^g are both scalar field dominated solutions with = 1. Comparing with 
point P4, Pg is the well-known scalar field dominated solution which exists for A^ < 6. 
TABLE 1 has shown that this scalar field dominated solution is a later-time attractor in 
the presence of a barotropic fluid if we have A^ < 87 and A=i,(i/* > 0. This solution will 
give an accelerating universe if A^ < 2 and X^:df^ > 0. For example, /(A) = ^ — 
corresponds to V^(0) = [cosh{l(i>)]" ■ Obviously we have A,, = ±|n(T| and df^ = The 
scalar field dominated solution with potential ^(0) = [cosh{lij})]" ^ late-time attractor if 
n^cr^ < 87 and > 0. In addition, this solution admits an accelerating expansion of 
universe if n^a^ < 2 and > 0. Noted that the point Pg means two stable critical points 

(Ac = ±\na\,Xc = ^^\na\,yc = ^Jl^^n?a^) in this case . 

Pio is the scaling solution where neither the scalar field nor the barotropic fluid entirely 
dominates the universe. Pio is a stable node for 247^/ (97 — 2) > A^ > 87 and X^df^ > 
and a stable spiral for A^ > 247^/(97 — 2) and X^:df^: > 0. So Pg and Pio can not be stable 
simultaneously. The scaling solution has drawn a lot of attentions since it can alleviate 
the coincidence problem of dark energy. Many potentials have been proposed to give a 
scaling evolution regime[85, 86, 52-65]. Here we give a sufficient condition for a potential 
to possess a scahng solution, that is, as long as /(A) equals zero for one or more values of 
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A(= A*) and these A* also satisfy Eq.(17), then there must exist a scaling solution with fi^ = 
37/A^. Obviously many potentials which satisfy this condition exist, such as the potential 
^(^) = [cosh(l,p)]" which corresponds to /(A) = ^ - the potential V^(0) = ^^^^-0.4.)^ which 
corresponds to /(A) = ^ + f and so on. Our condition includes the potentials in Ref[49] where 
the authors found that every positive and monotonous potential which was asymptotically 
exponential yielded a scaling solution. Our result is also not contradiction to the statement 
in literature [66, 67] where they assumed a scaling solution like Piq and found the potential 
was unique the exponential form. This exponential potential is explicitly figured out from the 
assumption and the evolution of universe with this potential is always the scaling solution(see 
Eq.(18) in literature [66]) while Piq being a stable point means that all the evolution of the 
universe with a class of potentials which satisfy Eq.(17) will all approach the scaling solution 
finally. It is just an asymptotic behavior at late time. Unfortunately, for the scaling solution 
of Pio, the state equation of dark energy equals Wm and therefore there does not exist the 
accelerating expansion if Wm is larger than zero. However, authors had obtained the exact 
quintessence potential V{^) = i^p^Jy^sm/i(^^^^^^^)]"2(i+«'<^)/Kn-«'^), which 

admited a scaling solution with 7^ Wm and VL^ 7^ 0[66]. With this potential, in principle, 
we can obtain a scaling solution with an accelerating expansion of the universe. 

Finally, let us consider the point P4, which is a de- Sitter-like dominant attractor with 
Vtif, = 1 and = —1. The condition for P4 being a stable point is that the value of /(A) 
when A = must be larger than zero(i.e., /(O) > 0, see appendix for details). So generally 
speaking, P4 and P9(or Pio) also can not be stable simultaneously. However, there may 
exist the possibility for some potentials that their values at A = is larger than zero but 
equals zero for some others A,,, (A* 7^ 0), then this region of A in the phase space of the three 
dynamical autonomous system will lie in the basin of the attractor Piq. That means, in this 
case, there can exist two stable critical points simultaneously, but this is not to say that 
the universe can evolve continuously from one stable critical point to another one. Based 
on this fact, the author proposed a scenario of universe which could evolve from a scaling 
attractor to a de-Sitter-like attractor by introducing a field whose value changed a certain 
amount in a short time[47]. In fact, we can also obtain these two asymptotical evolutions 
if the potential V{(j)) can be approximated to two different potentials when (j) evolves to 
different range, one admits the scaling solution and another one admits the de-Sitter-like 
solution[53, 58, 61]. For these potentials, the exit of the cosmological evolution from one 
attractor solution to another attractor is quite natural, but the explanation of why we have 
these special potentials is not quite natural. 
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5 Conclusion 



In this paper, we extend the autonomous dynamical system analysis of the canonical scalar 
field from 2-D to three-dimension by considering the potential parameter F as a function of 
another potential parameter A. There are ten critical points in all: three of these points 
(-^3,5,6) are general points which are possessed by all quintessence models regardless of the 
form of potentials and the rest points, with their existence or/and stability, are closely 
connected to the concrete potentials. We surprisingly find that, apart from the exponential 
potential, there are a large number of potentials which can give the scaling solution when 
the function /(A)(= r(A) — 1) equals zero for one or some values of A and the parameter A 
satisfies the condition Eq.(16) or Eq.(17) at the same time. We give the explicit expression 
to derive these potentials V{(f)) from /(A). We find that, if some conditions are satisfied, 
the de-Sitter-like dominant point P4 and the scaling point Pg{ or Pio) can simultaneously 
be stable, but Pg and Pio can not be stable at one time. As we have seen, the autonomous 
dynamical systems analysis is a very powerful tool which helps us extract useful cosmological 
information without solving the complicated background equations. Our method extends 
the analysis from two-dimensional autonomous dynamical system to three-dimension, which 
makes us be able to research a large number of potentials beyond the exponential potential. 
This method is quite effective and may be applied to a broad class of dark energy models 
studied in literature [46], including coupled quintessence, (coupled-) phantom scalar field, k- 
essence and even generalized background oc p^^. 

However, we should point out that our approach also has its drawbacks. First, as we have 
mentioned above: the approach can not be applied for the potentials for which the function 
r = VV" /{Vy can not be written as an explicit function of the variable A. Second, the 
variable A is undefined if the potentials vanish at its minimum, so the approach can not 
be applied for the potentials which vanish at its minimum. But, in despite of the second 
problem, it is actually not a fatal drawback. On one hand, the minimum of a potential 
is always associate with the late-time cosmological dynamics (future attractors). It is quite 
easy to discuss this special equilibrium point separately if we know a given potential has 
minimum(It is usually not difficult to find out the minimum of a given function). On the 
other hand, we can still use our approach to analyze all the critical points Pi_io since A is 
well-defined around these critical points. We take the potential V{(f)) = Vo[cos/i(a0) — 1] for 
example, this potential has a minimum value at = 0(A has no definition at = 0). We 
can investigate the critical point corresponding to this minimum separately. The explicit 
function about this potential is /(A) = |(^ — 1). Obviously, the point corresponding to the 
potential's minimum does not appear in the Table 1. We can still discuss the properties of 
the critical points Pi_io even if the variable A has no definition sometime. 
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Appendix 



In section 3, we pointed out that if the eigenvalues of Jacobi matrix had one or more eigenvalues 
with zero real parts while the rest of the eigenvalues had negative real parts, then linearization fails 
to determine the stability properties of this critical point. From TABLE 1 we realize that point P4 is 
just such point, so in this Appendix we will show you that how we get the stable condition of P4 from 
the center manifold theorem. The point P4 is (Ac = 0,Xc = 0,yc = 1) and its three eigenvalues 
are (0, —3, —3(1 + Wm))- Firstly, we transfer P4 to P4 (Ac = 0, Xc = 0, Fc = 1/c — 1 = 0) for 
convenience. In this case, Eqs.(5-7) can be rewritten as: 



dN 



-V6\^f{X)x (21) 



dx 
dN 
dY 



-3x + 



■Wm)--{l+W^)xY^-3il+Wr,^)xY (22) 



-3{l+w^)Y-^VQXxiY+l) + ^{l-w^)x''Y-^Y^-^{3+Wrn)Y^+^{l- 



-Wm. X 



(23) 



Noted that {X,x,Y} in Eqs. (21-23) are very small variables around point (Ac = 0,Xc 
0,Yc = 0. So the function /(A) in Eq.(21) should be taken the Taylor series in A: /(A) 
/(O) + p{0)X + ^A^ + where /"(O) is the value of when A = 0. 



., where /"(O) is the value of when A 

We can write down the Jacobi matrix A of dynamical system Eqs. (2 1-23): 



A 





rVe -s 








-3{l + w^) 
The eigenvalues of A and the corresponding eigenvectors are: 



{0, [1, ^,0]}; 



{-3, [0,1,0]}; {-3il + 3wm), [0,0,1]} 



(24) 



(25) 
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Let A4 he a matrix whose columns are the eigenvectors of A, then we can write down Ai 
and its inverse matrix T: 





1 










1 








M = 


6 


1 





T = M-' = 


6 


1 













1 










1 



(26) 



Using the similarity transformation T we can transform A into a block diagonal matrix, that 
is, 

^00 

T AT-^ =0-3 

-Sd+Wrr,) 



Al 

A2 



(27) 



where all eigenvalues of Ai have zero real parts and all eigenvalues of A2 have negative real 
parts. We put a change of variables: 



(28) 



" A' " 




' A " 




A 


x' 


= T 


X 




-fX + x 






Y 




Y 



Then we can rewrite the dynamical system Eqs. (21-23) in the form of new variables: 

dX' _ dX 
dN~ dN 

dX 



fiiX'^x'X) 



dx' 
dN 



dx 



6 dN dN 
dY' dY 



+ — = MX',x'X) 

h{x\x'X) 



(29) 

(30) 
(31) 



dN dN 

the detail forms of /i(A', x', F'), /2(A', x', F'), /3(A', x', Y') are easily obtained after we substi- 
tute the transformation A = A', x = ^A' + x' and Y = Y' into the right hand of Eqs. (21-23). 
According to the center manifold theorem, the stable condition of dynamical system Eq.(21- 
23), i.e., the stability of P4 will be finally determined by the following simple reduced system: 



-AV(0) 



(32) 



/(O) is the value of function /(A) at A = 0. This simple one-dimensional dynamical system 
Eq.(32) is stable if /(O) > 0. 

So we conclude that P4 is a stable de-Sitter-like dominant attractor when /(O) > 0, just 
as shown in TABLE 2. 
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